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A  series  of  high  power,  high  efficiency  Xa-band  and  W-band  gyroklystron  experiments  has  been 
conducted  recently  at  the  Naval  Research  Laboratory  (NRL).  Stagger  tuning  of  the  cavities  for 
bandwidth  enhancement  is  commonly  used  in  the  conventional  multicavity  klystrons.  The  desired 
stagger  tuning  is  usually  achieved  via  mechanical  tuning  of  the  individual  cavities.  However,  in  the 
millimeter  wave  regime,  particularly,  in  the  case  of  the  high  average  power  operation,  it  is  desirable 
to  be  able  to  achieve  the  required  stagger  tuning  by  design.  The  NRL  gyroklystron  experiments 
explored  the  tradeoffs  between  power,  bandwidth,  efficiency,  and  gain  to  study  the  effects  of  large 
stagger  tuning  in  millimeter  wave  without  resorting  to  mechanical  tuning  of  the  cavities.  Both, 
A^-band  and  W-band,  experiments  demonstrated  a  record  power-bandwidth  product.  The  success  of 
the  experiments  was  due  in  large  part  to  a  battery  of  improved  large-signal,  stability,  and  cold  test 
codes  employed  in  the  modeling  and  design  stage.  Theoretical  models  that  provide  the  basis  for 
these  design  simulation  tools  and  the  design  methodology  will  be  presented.  ©  1999  American 
Institute  of  Physics.  [S1070-664X(99)92905-1] 
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I.  INTRODUCTION 

Gyroklystrons  in  Ka- band  and  W-band  are  attractive 
sources  for  the  next  generation  of  millimeter  wave  radars.1 
This  is  because1  gyroklystrons  are  capable  of  operating  at 
both  high  peak  power  and  high  average  power  with  high 
efficiency  at  these  frequencies.2  The  instantaneous  band¬ 
width  of  gyroklystrons  is  determined  by  the  quality  factor 
(Q)  of  the  cavities,  and  especially  by  the  Q  of  the  output 
cavity.  This  leads  to  narrow  bandwidth,  since  the  Q  is  typi¬ 
cally  on  the  order  of  i50-300.  However,  it  is  possible  to 
increase  the  bandwidth  by  employing  a  multicavity  gy¬ 
roklystron  circuit,  with  each  cavity  tuned  to  a  different  fre¬ 
quency  (stagger  tuning).  Stagger  tuning  of  the  cavities  for 
bandwidth  enhancement  is  commonly  used  in  conventional 
multicavity  klystrons.  The  desired  stagger  tuning  is  usually 
achieved  via  mechanical  tuning  of  the  individual  cavities. 
However,  in  the  millimeter  wave  regime,  particularly,  in  the 
case  of  the  high  average  power  operation,  it  is  desirable  to  be 
able  to  achieve  the  required  stagger  tuning  by  design.  There¬ 
fore,  the  purpose  of  our  investigations  is  to  develop  and  vali¬ 
date  the  simulation  models,  the  design  methodology,  and  the 
manufacturing  procedure  to  achieve  this  goal  for  high  power 
gyroklystrons.  In  this  paper  we  describe  the  simulation  mod¬ 
els  and  the  design  procedure  to  address  the  relevant  issues. 

A  series  of  experiments  has  been  conducted  recently  at 


the  NRL  to  study  the  effects  of  large  stagger  tuning  on  the 
performance  of  the  gyroamplifiers.  These  experiments  ex¬ 
plored  the  tradeoffs  between  power,  bandwidth,  efficiency, 
and  gain  due  to  circuit  geometry  and  operating  parameter 
variation.  The  A^-band  experiments,  using  a  three-cavity3 
circuit  with  heavy  stagger  tuning,  demonstrated  more  than 
double  the  bandwidth  of  the  two-cavity  circuit,4  with  only  a 
modest  decrease  in  efficiency.  The  W-band  experiments 
achieved  a  record  power-bandwidth  product  in  a  four-cavity 
gyroklystron  circuit  with  large  stagger  tuning.5,6  For  all  ex¬ 
periments  the  results  are  in  excellent  agreement  with  simu¬ 
lations. 

Extensive  theoretical  work  has  produced  a  hierarchy  of 
nonlinear  models  to  describe  the  beam-wave  interaction  in 
the  gyroklystron  circuit.  The  validity  of  these  models  and 
their  accuracy  will  be  discussed  in  Sec.  II.  Ensuring  stable 
operation  is  a  central  step  in  the  gyroklystron  design  proce¬ 
dure.  Issues  relevant  to  self-oscillations  in  the  drift  sections 
of  the  gyroklystron  circuit  will  be  discussed  in  Sec.  HI.  The 
design  procedure  of  the  input  coupler  will  be  described  in 
Sec.  IV.  Finally,  in  Sec.  V  the  simulation  results  are  com¬ 
pared  with  the  experimental  data  and  the  methodology  de¬ 
veloped  for  designing  the  high  average  power  W-band  gy¬ 
roklystron  is  briefly  discussed. 
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II.  THEORETICAL  MODEL 

The  combination  of  the  phase  bunching  due  to  the  cy¬ 
clotron  maser  instability  in  a  cavity  and  the  ballistic  phase 
bunching  in  drift  tubes  makes  the  gyroklystron  capable  of 
operating  at  high  efficiency,  high  power,  and  high  gain.  Sev¬ 
eral  models  describing  the  gyroklystron  interaction  have 
been  developed  over  the  years,  ranging  from  an  analytical 
“point-gap”  approach7  to  a  multimode,  self-consistent  ap- 
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proach  implemented  in  a  simulation  code.8  A  complete  de¬ 
scription  of  the  variety  of  theoretical  models,  approaches, 
and  methods  is  beyond  the  scope  of  this  paper.  In  this  section 
we  will  describe  briefly  the  theoretical  models,  which  have 
been  implemented  in  the  large  signal  codes  used  at  the  Naval 
Research  Laboratory  (NRL)  for  designing  the  gyroklystron 
experiments. 

We  started  from  the  model  described  by  P.  Latham 
et  al9  The  assumptions  in  that  model  were  dictated  by  the 
intended  gyroklystron  application,  drivers  for  the  next  gen¬ 
eration  of  linear  colliders  for  particle  physics  research.  The 
collider  application  required  a  narrow  band,  low  average 
power  gyroklystron  capable  of  producing  a  very  high  peak 
power  with  high  efficiency.  This  was  achieved  by  using  a 
moderately  relativistic  500  kV  and  200  A  electron  beam.10 
The  NRL  gyroklystron  development  is  driven  by  a  radar  ap¬ 
plication,  which  requires  operation  at  high  average  power 
and  wide  bandwidth,  while  maintaining  relatively  high  peak 
power  and  efficiency.11  This  can  be  achieved  by  using  a 
weakly  relativistic  electron  beam  (50-70  kV  and  4-10  A)  in 
a  stagger-tuned  multicavity  circuit. 

It  turns  out  that  the  model  in  Ref.  9,  when  augmented 
with  additional  capabilities  required  for  the  design  of  multi¬ 
cavity  gyroklystron  with  stagger  tuning  for  bandwidth  en¬ 
hancement,  contains  essentially  all  of  the  effects  needed  for 
an  adequate  description  of  the  beam-wave  interaction  in  a 
gyroklystron.  Specifically,  a  more  accurate  model  of  the  in¬ 
put  coupler  and  a  model  for  calculating  the  dielectric  loading 
of  the  intermediate  cavities  have  been  added.  In  addition,  it 
turned  out  that  the  procedure  for  designing  a  gyroklystron 
amplifier  used  in  Ref.  9  is  not  suitable  for  effective  evalua¬ 
tion  of  the  bandwidth  of  the  circuit.  This  deficiency  was 


removed  as  well.  All  of  these  efforts  resulted  in  a  new  pro¬ 
cedure  for  designing  wide  band  gyroklystrons  and  to  an  im¬ 
proved  version  of  the  large  signal  code,  which  we  named 
MAGYKL.  Here  we  will  describe  only  the  additions  to  the 
model,  keeping  the  same  notation  as  Ref.  9. 

We  start  with  the  evolution  of  the  complex  amplitude  in 
each  cavity  in  MAGYKL,  which  is  described  by  a  time- 
dependent  equation 
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where  a  is  the  complex  normalized  amplitude  of  the  fields; 
Atu=2/{Re{o>c}~a)}/w  is  the  normalized  frequency  shift,  Qt 
is  the  load  quality  factor  of  the  cavity,  co  is  the  drive  fre¬ 
quency,  and  (oc  is  the  cold  resonant  frequency,  J  is  the  beam 
current  density,  WEM  is  the  field  stored  energy,  Ec  is  the  cold 
cavity  electric  field,  J  is  current  density,  T^lir/a)  is  the 
period,  and  r  is  the  slow  time12  normalized  to  Qi/o).  Note 
that  in  the  original  formulation  (Ref.  9),  Eq.  (1)  does  not 
contain  the  time  derivative.  We  will  come  to  this  point  later 
when  addressing  the  design  procedure. 

The  cavities  are  modeled  by  a  series  of  straight  uniform 
sections  with  abrupt  radial  discontinuities  at  the  junctions. 
The  fields  in  each  section,  expanded  as  a  radial  series  of  TE, 
TM,  and  TEM  modes,  are  determined  through  a  scattering 
matrix  solution.13  In  order  to  accurately  predict  the  amplifier 
bandwidth  of  the  amplifier,  the  input  cavity  model  used  in 
the  past  was  modified  to  include  a  frequency-dependent  elec¬ 
tric  field  amplitude  in  the  drive  cavity,  as  dictated  by  the 
resonant  frequency  and  total  Q  of  the  input  cavity.  The  field 
amplitude  in  the  drive  cavity  is  given  by 
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where  Pin  is  the  input  power,  which  is  assumed  to  be  fre¬ 
quency  independent,  G  external  is  the  external  quality  factor, 
Glosses  is  the  quality  factor  which  includes  ohmic  and  any 
other  losses  from  the  input  cavity,  such  as  field  leakage 
through  the  beam  tunnel,  vL  and  vz  are  the  perpendicular  and 
axial  electron  velocities,  £=zo)/c  is  the  normalized  distance, 
and  the  angle  brackets  denote  an  average  over  particles.  In 
the  derivation  of  Eq.  (2),  we  assumed  that  the  electric  field  in 
the  input  cavity  is  small  enough  that  the  beam -wave  interac¬ 
tion  can  be  evaluated  in  the  linear  regime. 

Before  describing  the  design  procedure,  a  few  remarks 
concerning  usage  of  “reduced  models”  are  in  order.  It  was 
originally  expected  that  the  model  used  in  Ref.  9  was  unnec¬ 
essarily  detailed  for  a  weakly  relativistic  and  low  current 
beam.  Therefore  a  reduced  model14  based  on  a  single  mode 
sinusoidal  or  Gaussian  axial  profile  for  the  “cold”  field 


could  be  used.  Although  reduced  models  of  gyroklystron  am¬ 
plifiers  are  still  useful  for  understanding  the  basic  operation 
of  these  devices,  they  are  not  adequate  as  a  design  tool  if  the 
combination  of  the  highest  possible  bandwidth  and  efficiency 
is  desired.  We  illustrate  the  limitations  of  the  reduced  models 
by  comparing  the  results  produced  by  these  models  with  our 
improved  large  signal  code,  MAGYKL.  We  will  use  the  four 
cavity  W-band  gyroklystron  with  stagger  tuning  as  an  ex¬ 
ample.  If  the  “cold  field”  in  the  cavities  is  represented  only 
by  one  radial  mode,  TE01,  and  the  axial  profile  is  sinusoidal, 
the  output  power  versus  input  power  (the  transfer  function) 
will  be  significantly  different  from  the  transfer  function  pre¬ 
dicted  by  MAGYKL.  This  is  shown  in  Fig.  1 .  Similarly,  the 
power  versus  frequency  dependence  (bandwidth  curve)  at 
saturation  as  predicted  by  magykl  is  quite  different  from  the 
prediction  of  a  model  based  on  a  single  mode  sinusoidal  or 
Gaussian  axial  profile  for  the  “cold”  field. 
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FIG.  1.  Efficiency  versus  input  power  in  the  case  of  four-cavity  gy- 
roklystron.  The  four-cavity  structure  with  the  following  parameters:  Cavity 
1  fc= 93.4064  GHz,  Gext=150,  £ioss=800;  Cavity  2  fc  =  94.2063  GHz, 
2ioaded=175;  Cavity  3  fc= 93 .2836  GHz,  Q =  175;  Cavity  4  fc 
=  93.8874  GHz,  (2 loaded =  162.43;  electron  beam  parameters  are:  y 
=  65kV,  I~6A  a  =1.5,  guiding  center  radius  rg  =  0.081  cm,  zero  velocity 
spread.  Magnetic  field  was  optimized  for  maximum  efficiency,  B 
=  36.8kG.  , 


A.  Design  procedure 

The  process  of  using  magykl  for  designing  a  gy- 
roklystron  is  a  multistep  procedure,  which  includes  the 
choice  of  an  electrodynamical  structure,  the  electron  beam 
and  drive  signal  parameters,  and  the  optimization  of  device 
performance  by  varying  the  magnetic  field.  The  determina¬ 
tion  of  the  electrodynamical  structure  consists  of  several 
steps.  First,  the  resonance  frequency,  g,  and  length  for  each 
cavity  are  chosen.  Based  on  this  information,  CASCADE13  is 
used  to  determine  approximate  cavity  shapes  and  sizes  and 
the  drift  tube  radii  that  yield  the  desired  frequencies  and  Q  s 
as  well  as  the  electric  field  profile,  Ec(r,z).  The  field  axial 
profile  and  the  amplitudes  and  phases  of  the  modes  together 
with  the  “cold”  frequencies  and  Q  s  are  used  in  the  large 
signal  code  magykl,  in  which  the  particle  equations  and  the 
wave  amplitude  Eq.  (1)  are  solved.  At  this  point  the  input 
cavity  parameters,  such  as  drive  frequency,  Pin,  gext,  and 
gloss,  the  beam  parameters,  such  as  current,  voltage, 
a~v±fvz ,  root-mean-square  (rms)  velocity  spread  and  the 
guiding  center  radius,  and  the  magnetic  field  value  are  cho¬ 
sen. 

The  amplitude  in  the  input  cavity  is  determined  by  Eq. 
(2)  via  iterations.  Usually,  three  to  five  iterations  are  suffi¬ 
cient.  The  number  of  the  time  steps  needed  for  the  ampli¬ 
tude,  a,  for  the  cavities  other  than  the  input  to  reach  a  steady- 
state  solution  is  usually  on  the  order  of  100.  The  overall 
efficiency  is  evaluated  using  this  steady-state  amplitude. 
Note  that  the  steady-state  solution  is  consistent  with  the  pre¬ 
scribed  total  quality  factor. 

In  the  original  version  of  the  large  signal  code,9  without 
the  time  derivative  term,  the  values  Q  and  the  magnitude  and 
phase  of  the  complex  amplitude  a  and  its  phase  and  Q  were 
used  to  determine  from  the  energy  balance  equation,  as 
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in  Eq.  (7)  in  Ref.  9.  This  approach  proved  to  be  computa¬ 
tionally  effective  for  the  cases  studied  in  Ref.  9.  However,  it 
turned  out  to  be  not  practical  when  stagger  tuning  was  em¬ 
ployed  for  both  efficiency  and  the  bandwidth  optimization. 
Adding  the  time  derivative  did  not  change  the  nature  of  the 
final  solution,  since  the  complex  amplitude  has  to  reach  a 
steady  state  in  order  to  satisfy  the  power  balance  equation. 
However,  the  numerical  implementation  of  solving  Eq.  (2)  is 
such  that  the  amplitude  and  phase  of  the  wave  are  adjusted 
from  one  time  step  to  another  automatically  until  the  power 
balance  equation  is  satisfied.  The  computation  is  repeated  for 
different  values  of  the  magnetic  field  in  order  to  optimize  the 
performance.  Note  that  the  drift  length  is  mainly  determined 
by  the  gain  or  efficiency  optimization. 

The  choice  of  the  guiding  center  radius  involves  addi¬ 
tional  considerations.  On  the  one  hand,  it  is  desired  to  make 
the  field  of  the  operating  mode  evanescently  small  in  the 
drift  region,  hence  a  small  drift  tube  radius  is  advantageous. 
On  the  other  hand,  the  beam  size  has  to  be  adequate  to  en¬ 
sure  maximum  coupling  to  the  electric  field  and  minimum 
interception  by  the  drift  walls.  Usually  a  distance  on  the 
order  of  three  Larmor  radii  from  the  drift  wall  is  considered 
to  be  safe.  Another  interesting  point  worth  mentioning  con¬ 
cerns  the  adequacy  of  the  small  Larmor  radius  approxima¬ 
tion  that  is  typically  used  for  describing  the  beam-wave  in¬ 
teraction  in  the  weakly  relativistic  regime.  We  find  that  if  the 
small  Larmor  radius  approximation  is  used,  the  frequency 
pulling  due  to  the  beam  loading  is  about  50  MHz  different 
from  magykl  calculations  (no  small  Larmor  radius  approxi¬ 
mation  is  used  in  Bessel  function  evaluation  in  the  equations 
used  in  magykl).  By  itself  it  is  a  tiny  error,  blit  considering 
the  fact  that  the  stagger  tuning  of  the  cavities  is  on  the  order 
of  100-150  MHz,  the  50  MHz  error  is  quite  significant. 

B.  Effects  of  the  dielectric  insert 

The  total  quality  factors  for  the  intermediate  cavities  are 
assumed  known  in  our  model.  In  practice,  in  order  to  realize 
the  required  cavity  losses,  a  ring  of  a  lossy  dielectric  is  in¬ 
serted  into  the  cavities,  magykl  requires  prior  calculation  of 
the  field  eigenmode  profiles  for  each  cavity,  which  are  nor¬ 
mally  obtained  from  the  scattering  matrix  code  CASCADE, 
which  does  not  allow  the  ceramic  insert  to  be  included  in  the 
computation.  Therefore  some  approximation  was  necessary 
to  model  the  actual  device.  Using  MAGYKL,  a  design  begins 
with  a  set  of  frequency  and  Q  values  for  each  cavity,  and  a 
known  dielectric  constant  for  the  ceramic  inserts. 

An  unpublished  code,  DIELEC,  developed  by  P.  Latham 
and  modified  by  B.  Levush  and  A.  Berman,  is  used  for  each 
cavity  to  find  dimensions  (radius,  length  of  the  dielectric 
ring),  dielec  is  a  more  advanced  scattering  matrix  code  that 
permits  (lossy)  ceramic  rings  in  the  description  of  the  cavity 
geometry.  Loss  is  included  through  the  use  of  complex  val¬ 
ues  of  the  dielectric  constant.  The  cavity  frequency  and  Q 
obtained  are  passed  as  inputs  to  MAGYKL.  Note  that,  when 
DIELEC  is  used,  the  number  of  modes  required  to  converge  on 
the  eigenfrequency  is  about  10.  However,  since  cascade 
must  be  used  for  the  field  profiles  passed  to  magykl,  ap¬ 
proximate  profiles  are  computed  with  the  assumption  that  the 
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cavity  electromagnetic  fields  are  excluded  from  the  dielectric 
ring  as  if  it  were  metal.  This  assumption  is  valid  if  the  real 
part  of  the  dielectric  constant  of  the  ceramic  is  large.  Cavity 
dimensions  for  this  equivalent  cavity  are  kept  the  same  as 
those  of  the  true  cavity,  except  that  the  radius  used  in  CAS¬ 
CADE  is  adjusted  slightly  to  match  the  required  resonant  fre¬ 
quency.  We  have  recently  adapted  magykl  to  take  fields 
from  dielec.  Detailed  description  of  the  model  will  be  pre¬ 
sented  in  other  publications.  The  adapted  code,  magykl-C, 
follows  a  similar  algorithm,  but  instead  of  requiring  an  ap¬ 
proximate  equivalent  cavity  to  be  constructed,  the  field  pro¬ 
file  is  taken  directly  from  DIELEC.  The  cavity  frequency  and 
Q  are  therefore  fully  consistent  with  the  cavity  and  dielectric 
dimensions,  so  that  a  level  of  approximation  has  been  re¬ 
moved  from  the  calculation. 

The  original  code  magykl  has  been  used  successfully  in 
the  design  of  NRL  gyroklystron  experiments  using  the  cavity 
field  approximation  described,  and  therefore  for  the  param¬ 
eter  regime  of  those  experiments  the  approximation  of  field 
exclusion  from  the  dielectrics  seems  most  likely  valid.  Here 
we  explore  this  validity  by  comparing  results  based  on  the 
original  model  with  the  more  complete  model,  and  extend 
the  comparison  to  other  parameter  regimes  of  experimental 
interest.  We  will  consider  two  cases  of  dielectric  properties: 
(A)  large  dielectric  constant,  high  loss  (e=  12.24-/3.67) 
and  (B)  large  dielectric  constant,  moderate  loss  (e=10.0 
—  /0.88).  MAGYKL  has  been  used  successfully  in  the  design 
of  devices  in  the  case  “A”  regime,  taking  the  field  profile  as 
if  the  ceramic  were  metal.  Case  “B”  requires  a  greater 
length  of  ceramic  insert,  necessary  to  achieve  the  same  over¬ 
all  absorption  to  achieve  the  design  Q.  The  DIELEC  field  pro¬ 
file  (“true”  cavity)  has  a  more  complex  structure,  and  the 
cascade  field  profile  used  by  magykl  is  not  as  good  an 
approximation  in  this  parameter  regime.  Figure  2  shows  the 
profile  of  the  cavity  field  for  the  radial  position  at  which  the 
electron  beam  travels,  for  both  the  true  cavity  and  the 
equivalent  cavity  for  case  B.  The  effect  of  the  extended  di¬ 
electric  region  is  apparent,  causing  a  plateau  in  the  field  pro¬ 
file.  This  deviation  from  the  equivalent  cavity  profile  is  suf¬ 
ficient  to  influence  the  electron  beam.  Figure  3  shows  the 
transfer  function  curves  calculated  by  both  magykl  and 
magykl-C  for  this  case.  It  can  be  seen  that  magykl  and 
MAGYKL-C  differ  distinctly  in  the  peak  efficiency.  Similarly, 
the  efficiency  versus  frequency  dependence  as  predicted  by 
magykl  and  magykl-c  is  different.  For  this  design,  it 
would  be  necessary  to  use  the  new  version  of  the  code  to  be 
confident  of  the  results.  However,  in  the  cases,  such  as  case 
A,  the  benefits  of  the  computation  speed  justify  the  proce¬ 
dure  of  using  th^  “cold”  field  axial  profile  in  MAGYKL  with¬ 
out  taking  into  account  the  effect  of  the  lossy  dielectric. 

III.  DRIFT  TUBE  DESIGN 

Gyroklystron  drift  spaces  ideally  provide  a  field-free  en¬ 
vironment  in  which  ballistic  bunching  of  the  electron  beam 
can  occur.  As  discussed  in  preceding  sections,  the  drift  tube 
radius  is  made  small  enough  to  ensure  that  the  cavity  oper¬ 
ating  mode  is  cut  off  inside  the  drift  tube.  However,  this  does 
not  guarantee  that  the  drift  tube  is  free  from  the  fields  asso- 
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FIG.  2.  The  cavities’  parameters  are  as  follows:  Cavity  1  fc 
-92.486  GHz,  <2ext=300,  (2 loss =  800;  the  drift  length  between  the  first  and 
second  cavity  0.64  cm;  Cavity  2  length  L2  =  0.5cm,  fc 

=  93.494  GHz,  f2ioaded=  170;  the  drift  length  between  the  second  and  third 
cavity  Lddft,2=  0.82  cm;  Cavity  3  length  L3  =  0.5cm,  /c  =  92.917  GHz, 
<2  loaded =  170;  the  drift  length  between  the  third  and  forth  cavity  3 
=  0.82  cm;  Cavity  4  fc= 93.3 15  GHz,  (2ioaded=338;  the  dielectric  permittiv¬ 
ity,  e  =  10.00—  z'0.88  has  a  slightly  lower  real  part,  permitting  greater  field 
penetration,  and  a  smaller  loss  tangent,  so  that  a  larger  region  of  dielectric  is 
necessary  to  achieve  the  same  value  of  design  Q.  This  was  achieved  with  the 
dielectric  inserts  of  width  Ldiel  cavl  =  Ldiel  cav2=0.18  cm.  The  solid  line  in  the 
figure  is  the  electric  field  profile  in  die  second  cavity  determined  by  cas¬ 
cade.  The  dashed  line  in  the  figure  is  the  electric  field  profile  in  the  second 
cavity  determined  by  dielec. 

dated  with  self-oscillation.  A  simple  metal  drift  tujbe, 
bounded  by  transitions  to  the  gyroklystron  cavities  on  either 
end,  can  function  as  an  unintentional  resonant  cavity  to 
modes  different  from  the  normal  operating  mode.  This  is 
particularly  true  with  gyroklystrons  operating  in  the  TE0i 
mode,  since  the  lower  order  TEn  mode  can  propagate  at  drift 


Pto(W) 


FIG.  3.  The  transfer  curve  predicted  by  magykl  (dielectric  insert  assumed 
not  perturbing  the  fields  in  the  cavity,  hence  cascade  fields  are  used)  and  by 
magykl-c  (dielec  fields  are  used).  The  electromagnetic  structure  parameters 
are  the  same  as  in  Fig.  2.  The  beam  parameters  are:  V=55kV,  7=4A,  a 
=  1.5,  rms  (Ayx  fvx)-  10%.  Magnetic  field  is  35.8  kG  and  the  drive  fre¬ 
quency  is  93.45  GHz. 
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tube  radii  and  frequencies  considerably  below  the  limiting 
values  for  the  TE01  mode.  These  types  of  instabilities  asso¬ 
ciated  with  inadvertent  drift  tube  resonators,  as  well  as 
backward-wave  gyro-instabilitips  in  the  drift  tubes,  become 
increasingly  troublesome  as  one  pushes  toward  higher  beam 
power. 

To  prevent  drift  tube  oscillations,  lossy  dielectric  mate¬ 
rials  (absorbers)  are  placed  inside  the  drift  spaces.15  The  ab¬ 
sorber  consists  of  a  continuous  liner  of  lossy  dielectric  that 
covers  most  of  the  drift  tube  length.  The  extreme  ends  of  the 
drift  tube  adjacent  to  the  cavities  are  metal,  to  prevent  power 
deposition  into  the  dielectric  from  leakage  fields.  Often  the 
dielectric  radius  is  made  somewhat  larger  than  the  surround¬ 
ing  metal,  to  help  prevent  charge  buildup  from  stray  elec¬ 
trons  from  the  beam.  The  procedure  for  analyzing  the  stabil¬ 
ity  of  a  candidate  drift  tube  consists  of  a  number  of  steps  that 
are  applicable  to  any  of  the  configurations,  although  the  de¬ 
tails  will  depend  on  the  specific  structure  under  investigation. 
Essentially,  one  first  determines  the  frequencies  where  oscil¬ 
lation  could  occur.  This  is  followed  by  an  evaluation  of  the 
attenuation  and  quality  factor  of  the  lossy  structure  at  these 
frequencies,  in  comparison  to  computed  start  oscillation  val¬ 
ues.  Finally,  if  needed,  the  physical  dimensions  of  the  drift 
tube  structure  can  be  adjusted  to  optimize  performance  in  an 
iterative  manner. 

To  illustrate  the  procedure,  consider  placing  a  liner  with 
a  radial  thickness  of  0.8  mm  and  a  complex  dielectric  con¬ 
stant  of  e  =  12.24-  Z3.67  around  the  1 .75  mm  drift  tube.  Pro¬ 
vided  that  the  length  of  the  dielectric  is  at  least  several  wave¬ 
lengths  long,  one  can  ignore  edge  effects  and  model  the 
structure  as  being  infinitely  long. 

Cyclotron  maser  instabilities  are  most  likely  to  occur  at  a 
frequency  where  the  dispersion  curve  (a  plot  of  frequency 
versus  axial  wave  number)  of  the  drift  structure  intersects 
with  the  dispersion  curve  for  Doppler  upshifted  cyclotron 
waves  supported  by  the  electron  beam.  When  a  dielectric 
liner  is  placed  into  the  drift  tube,  two  effects  occur.  First,  for 
modes  that  can  propagate  in  the  empty  guide  at  the  fre¬ 
quency  of  interest,  the  presence  of  losses  creates  a  complex 
propagation  factor  in  the  axial  direction  /3  =  kz-ikzi ,  where 
the  axial  variation  of  fields  is  assumed  to  follow  the  func¬ 
tional  form  e~jPz.  The  term  kzi  represents  an  attenuation 
caused  by  the  losses.  The  second  effect  of  the  dielectric  is  a 
significant  change  in  kz ;  modes  that  were  originally  heavily 
cut  off  at  the  frequency  of  interest  (kzi> 0,  kz= 0)  can  often 
propagate  as  a  result  of  the  addition  of  dielectric.  For  the 
lossy  dielectric  liner  to  be  effective,  it  is  important  that  these 
newly  liberated  modes  receive  sufficient  attenuation  as  well. 
In  the  same  manner,  the  addition  of  dielectric  also  modifies 
the  functional  form  of  the  dispersion  curve  and  changes  the 
instability  frequencies.  The  modified  dispersion  curves  for 
what  originally  were  the  TEn,  TMU,  and,TE12  modes  are 
shown  in  Fig.  4,  and  are  now  indicated  by  HEn,  HMn,  and 
HE12.  This  latter  notation  indicates  that  the  modes  are  now 
hybrid  (both  Ez  and  Hz  components  exist).  The  new  HEn 
backward-wave  intersection  frequency  with  the  beamline  is 
now  at  51.8  GHz,  while  the  forward-wave  interaction  with 
this  mode  is  eliminated.  However,  the  HMn  and  HE12  curves 
have  been  pulled  down  to  have  intersections  at  56.0  GHz  for 
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FIG.  4.  The  dispersion  curves  are  plotted  for  the  TEU>  TMn,  and  TE12 
modes  for  a  =1.75  mm.  The  dispersion  curve  for  cyclotron  waves  on  a 
relativistic  beam  is  given1  by  a)=Q,+kz vz,  where  f 1  is  the  relativistic  cy¬ 
clotron  frequency,  kz  is  the  axial  wave  number  (wave  propagation  factor), 
and  vz  is  the  axial  beam  velocity.  It  is  represented  by  the  straight  line  in  the 
figure,  for  the  case  of  a  65  kV  beam  with  a  36.8  kG  magnetic  field  and  a 
=  1.5.  The  modified  dispersion  curves  for  what  originally  were  the  TEU, 
TMn,  and  TE12  modes  are  also  shown  in  the  figure,  and  are  now  indicated 
by  HEn,  HMU,  and  HE12.  This  latter  notation  indicates  that  the  modes  are 
now  hybrid  (both  Ez  and  H.  components  exist). 


the  backward  wave  HMn,  75.2  GHz  for  the  backward-wave 
HE12,  and  119.2  GHz  for  the  forward- wave  HE12.  Having 
identified  the  relevant  frequencies  for  beam-wave  intersec¬ 
tions  (and  thus  potential  instabilities),  one  can  now  evaluate 
the  stability  of  the  loaded  structure. 

To  better  quantify  the  effect  of  dielectric  on  each  of  the 
modes,  one  must  realize  that  is  often  not  the  attenuation  per 
unit  length  that  is  critical  for  suppressing  oscillation,  but 
rather  the  quality  factor  of  the  inadvertent  cavity  that  can 
exist  within  the  drift  tube.  This  quality  factor  will  determine 
whether  a  particular  beam  power  will  trigger  oscillation 
since  stability  scales  as  the  product  of  beam  power  with  Q.  If 
the  transitions  on  both  ends  of  the  drift  tube  produced  total 
reflection  (a  worst-case  situation),  the  dielectric-induced 
quality  factor  of  such  a  cavity  would  be  Q  =  Qn=kzlkzi .  In 
Fig.  5  the  quality  factors  of  the  modes  are  plotted  versus 
frequency.  The  potential  instability  frequencies  from  Fig.  4 
are  indicated  with  the  circles.  One  can  see  that  quality  factors 
at  the  HEn,  HMU,  and  the  lower  HE12  instability  frequen¬ 
cies  are  below  7,  which  indicates  extremely  effective  load¬ 
ing.  Since  these  values  of  Qn  are  massively  below  the  empty 
beam  tunnel  starting  Q  for  the  design  current  (about  200), 
drift  tube  stability  at  these  intersection  frequencies  is  guar¬ 
anteed.  Because  the  gyroklystron  amplifier  will  be  used  in 
radar,  it  must  be  unconditionally  stable  at  the  operating  point 
in  the  cavities  as  well  as  the  drift  tubes.  Thus,  a  linear  analy¬ 
sis  was  used  to  determine  stability  conditions  in  each  part  of 
the  circuit.16  More  detailed  stability  analysis  is  out  of  the 
scope  of  this  paper  and  will  be  published  elsewhere  in  the 
future. 
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FIG.  5.  The  quality  factors,  Q  =  kzlkzi,  of  the  modes  are  plotted  versus 
frequency.  The  potential  instability  frequencies  from  Fig.  4  are  indicated 
with  the  circles. 


IV.  DESIGN  OF  THE  INPUT  COUPLER 

The  input  coupler  has  to  satisfy  the  following  require¬ 
ments.  First,  the  bandwidth  of  the  coupler  should  be  enough 
to  meet  the  demands  of  the  gyroklystron  bandwidth.  Second, 
the  operating  mode  should  be  TE01  and  there  should  be  high 
degree  of  mode  purity.  Third,  the  amount  of  power  coupled 
into  the  TE01  mode  should  be  maximized.  To  meet  these 
requirements  we  choose  a  coaxial  slotted  structure  for  the 
input  coupler,17  shown  in  cross  section  in  Fig.  6.  The  coupler 
has  the  advantage  of  a  single  guide  feed,  good  cavity  mode 
purity,  minimal  ohmic  loading  of  the  cavity  walls,  and  it  can 
be  designed  with  a  low  loaded  quality  factor. 

The  input  wave  guide  excites  the  TE411  mode  in  the 
outer  cavity.  Power  is  coupled  through  four  axial  slots  to  the 
inner  cavity  TE011  mode.  We  employed  in  the  design  effort 
two  general  three-dimensional  (3D)  codes,  HFSS18  and 


input 

power 


FIG.  6.  hfss  simulation  of  coaxial  drive  cavity  in  cross  sectional  view 
showing  |£|2. 


ARGUS.19  These  codes  are  the  only  practical  tool  to  investi¬ 
gate  the  effects  of  cylinder/coax  centerline  offsets,  coupling 
slot  axial  offsets,  nonequal  slot  azimuthal  spacing,  slot  rota¬ 
tions,  cavity  wall  distortions,  and  other  nonideal  configura¬ 
tions.  These  codes  are  viewed  as  complementary  in  their 
characteristics;  each  with  particular  advantages  that  can  be 
exploited.  The  tetrahedral  modeling  elements  employed  by 
HFSS  conform  nicely  to  structures  with  curved  surfaces,  al¬ 
though  small  “facets”  must  be  used  for  high  accuracy.  This 
results  in  a  large  number  of  unknowns,  leading  to  large  prob¬ 
lems  that  strain  computational  resources.  Geometric  param¬ 
eterization  is  also  difficult  due  to  the  adaptive  meshing,  since 
a  geometry  change  results  in  a  different  initial  mesh.  On  the 
other  hand,  the  ability  to  transmit  a  wave  through  the  struc¬ 
ture  at  a  particular  frequency  closely  mimics  typical  testing 
or  operation  scenarios.  ARGUS  clearly  identifies  the  cavity 
eigenmodes,  in  contrast  to  HFSS,  where  decomposition  of  the 
electric  field  is  required  to  determine  mode  content.  The 
structured  grid  of  ARGUS  results  in  a  “stairstep”  approxima¬ 
tion  of  curved  surfaces,  requiring  a  dense  grid  for  adequate 
resolution.  ARGUS  is  somewhat  faster*  since  it  can  be  run  on 
Cray  supercomputers,  although  the  overall  solution  time  is 
comparable  to  that  using  hfss,  since  multiple  runs  are  re¬ 
quired  for  Q  determination.  Large  problems  can  be  handled 
with  ARGUS,  albeit  with  long  solution  times.  However,  pa¬ 
rameterization  is  relatively  easy  with  a  structured  grid.  Fi¬ 
nally,  the  real  utility  of  using  both  codes  for  coupler  model¬ 
ing  is  confidence  in  the  results.  The  same  answers  are 
obtained  using  two  completely  different  solution  approaches. 
Both  codes  were  used  with  various  methods  of  solution  to 
model  many  issues  concerning  the  input  coupler,  including 
resonant  frequencies,  modes,  g’s,  resistive  loading,  mode 
purity,  parameterization,  and  manufacturing  tolerances.-  In 
particular,  ARGUS  and  hfss  were  used  to  model  the  Q  and 
resonant  frequency  for  the  operating  mode  of  the  input  cou¬ 
pler.  Excellent  agreement  (0.1%)  was  achieved  between  the 
two  codes. 

Several  techniques  were  employed  with  HFSS  to  deter¬ 
mine  the  resonant  frequency  and  Q  of  the  TE01  mode  in  the 
cylindrical  cavity.  The  cavity  resonant  frequency  and  Q  were 
determined  from  the  frequency  dependence  of  stored  energy. 
The  resonant  frequency,  fr ,  is  the  frequency  at  maximum 
stored  energy  and  Q=frf(fu~fi)*  where  fu  and/j  are  the 
upper  and  lower  frequencies  where  the  stored  energy  is  half 
of  the  maximum  value.  Specifically,  E*E  was  calculated 
over  the  cavity  volume,  using  the  “fast-frequency  sweep” 
option  which  uses  an  asymptotic  wave  form  evaluation  to 
calculate  a  “bandwidth”  solution  from  a  full  field  solution  at 
a  single  frequency.  Excellent  agreement  with  ARGUS  results 
was  obtained  with  a  sufficiently  high  number  of  polygon 
segments  approximating  circular  boundaries,  ~2°/facet.  A 
second  method  “de-embeds”  Sn  at  the  drive  port  to  find  the 
“detuned  short”  position.  The  resonant  frequency  and  Q  are 
calculated  from  the  phase  of  the  reflection  coefficient.  HFSS 
has  also  been  used  to  model  a  ‘  ‘transmission’  ’  measurement 
by  launching  the  TE01  mode  in  the  cylindrical  cavity  by  eva¬ 
nescent  coupling  through  the  drift  tube  and  “receiving”  at 
the  rectangular  port. 

ARGUS  is  used  in  the  eigenvalue  solver  mode.  In  prin- 
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ciple,  all  relevant  inodes  and  frequencies  can  be  determined. 
Since  the  entire  problem  volume,  including  conductors,  is 
meshed,  a  sector  of  the  configuration  is  typically  modeled  for 
reasonable  problem  size  and  solution  time.  Resonant  fre¬ 
quencies  and  Q  s  in  ARGUS  are  determined  from  a  variation 
of  the  “detuned  short”  method  due  to  Kroll  and  Yu20  requir¬ 
ing  multiple  runs  as  the  length  of  the  input  wave  guide  is 
varied.  The  method  of  geometry  input  in  ARGUS  is  more 
conducive  to  parametric  variation  than  for  HFSS,  once  the 
basic  geometry  has  been  set.  HFSS,  on  the  other  hand,  can  be 
more  quickly  set  up  when  major  geometry  changes  are  in¬ 
volved.  Both  HFSS  and  ARGUS  were  used  to  estimate  the 
mode  purity  of  the  TE011  circular  cavity  mode.  Performance 
of  the  input  coupler  is  found  using  both  HFSS  and  ARGUS  to 
be  quite  sensitive  to  dimensions,  and  the  dimensional  sensi¬ 
tivity  approaches  the  limits  of  ordinary  manufacturing  toler¬ 
ances.  Note  that  a  0.0003  cm  change  in  cylindrical  cavity 
radius  shifts  the  resonant  frequency  by  —100' MHz.  It  is 
important  that  the  slot  angular  orientation  be  within  about  1° 
of  the  normal  value  of  45°  with  respect  to  the  input  wave 
guide.  A  1°  slot  rotation  increases /by  —20  MHz  and  re¬ 
duces  Q  by  3%.  A  5°  slot  rotation  increases /by  —36  MHz 
and  reduces  Q  by  30%.  A  frequency  and  Q  within  ±100 
MHz  and  ±50,  respectively,  of  the  design  values  are  accept¬ 
able  from  overall  circuit  performance  considerations. 

V.  SUMMARY 

The  design  of  both  two  and  three  cavity  Ka- band  gy- 
roklystron  circuits  have  been  validated  experimentally.3,4 
The  design  of  the  W-band  gyroklystron  circuits  have  been 
validated  experimentally  with  record  performance  in  effi¬ 
ciency,  peak  power,  and  bandwidth,5,6  thereby  demonstrating 
the  validity  of  the  design  codes  and  the  design  methodology. 
For  example,  in  Fig.  7  (Fig.  8  from  Ref.  6)  the  measured 
efficiency  and  peak  output  power  versus  frequency  are  plot¬ 
ted  together  with  the  predicted  values.  The  theoretical  pre¬ 
dictions  are  in  very  good  agreement  with  the  measured  data. 
In  this  particular  example  the  W-band  gyroklystron  circuit 
performance  was  designed  to  mimic  the  electrical  perfor¬ 
mance  of  the  high  average  power  version  of  the  gy¬ 
roklystron.  Radar  applications  usually  require  devices  having 
high  peak  and  average  power  capabilities.  The  procedure  of 
the  circuit,  drift  section,  and  input  coupler  design  described 
above  has  been  used  in  designing  an  industrial  version  of  the 
NRL  W-band  gyroklystron,  namely  a  four-cavity  gy¬ 
roklystron  capable  of  operating  at  100  kW  peak  power,  with 
efficiency  more  than  20%  and  instantaneous  bandwidth  of 
600  MHz  and  capable  of  delivering  10  kW  average  power.11 
The  3D  codes  HFSS  and  argus  have  been  also  used  to  iden¬ 
tify  the  geometry  of  the  intermediate  cavities  that  provides 
the  assumed  eigenfrequencies,  quality  factors,  and  mode  pu¬ 
rity,  while  being  consistent  with  thermal  management  and 
material  constraints.  Furthermore,  magykl  was  used  to 
study  variation  in  circuit  performance  with  changes  in  circuit 
dimensions,  thus  establishing  manufacturing  tolerances.21 
The  thermal  management  issues  for  the  cavities  have  been 
addressed  by  using  the  3D  thermomechanical  code  ANSYS.22 
The  performance  of  the  high  average  power  W-band  gy- 
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FIG.  7.  Theoretical  (solid  line)  and  experimental  (black  circles)  peak  power 
and  efficiency  as  a  function  of  drive  frequency.  Cavity  l  fc= 93.00  GHz, 
Gex *=150,  2loss=800;  Cavity  2  /c  =  93.74  GHz,  eIoaded=108;  Cavity  3  fc 
=  93.14  GHz,  <2loaded=  107;  Cavity  4  fc  =  93.49  GHz,  0loaded=145;  Elec¬ 
tron  beam  parameters  are:  V=57.72kV,  7=4.2A,  a=\.5,  perpendicular 
velocity  spread,  AvL  lvy  =  9%.  Magnetic  field  was  optimized  for  maximum 
efficiency  and  bandwidth  above  600  MHz,  B  =  36.0kG.  The  3  dB  band¬ 
width  is  640  MHz. 

roklystron  will  demonstrate  the  extent  of  the  validity  of  the 
entire  design  methodology,  which  includes  design  tools  for 
both  electrical  and  thermomechanical  analysis. 
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